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Abstract 

We study a class of self similar processes with stationary increments belonging to higher order 
Wiener chaoses which are similar to Hermite processes. We obtain an almost sure wavelet-like 
expansion of these processes converging uniformly on every compact. This allows us to compute the 
pointwise and local Holder regularity of sample paths and to analyse their behaviour at infinity. We 
also provide some results on the HausdorfF dimension of the range and graphs of multidimensional 
0^ ■ anisotropic self similar processes with stationary increments defined by multiple Wiener integrals. 

AMS classification: 60G18, 60G17, 60G22, 42C40, 28 A 78, 28 A 80. 

p I . Key words: Self-similar processes, multiple Wiener-Ito integral, wavelet expansion. Holder regularity, 

(-H ! Hausdorff dimension. 
> 

S ■ 1 Introduction and background 

■ Self similar processes with stationary increments (SSSI processes), i.e. processes X which satisfy: 

m . vc> 0, x,t = c^Xt, 

id) 

y{s,t)eR^ Xt-Xs = Xt^s 
^; ^'^^ 

, for some positive H, have been studied for a long time due to their importance both in theory and 
^ I in practice. Such processes appear as limits in various renormalisation procedures [151 IMl [2H]- In 
addition, they are the only possible tangent processes [12] • In applications, they occur in various 
• I— I . fields such as hydrology, biomedicine and image processing. The simplest SSSI processes are simply 
^ [ Brownian motion and, more generally. Levy stable motions. Apart from these cases, the best known 
' such process is probably fractional Brownian motion (fBm), which was introduced in |14] and popu- 
larized by |16j . A construction of SSSI processes that generalizes fBm to higher order Wiener chaoses 
was proposed in [T9]. These processes read 



VtGM+ Xt= [ htixi,...,Xd)dB^,...dB^^ 



where ht verifies: 

1. ht £ 

2. Vc> 0, hct{cxi,...,cxd) = c^'iht{xi,...,Xd), 

3. Vp > 0, ht+p{xi, ...,Xd) - ht{xi, ...,Xd) = hp{xi - t, ...,Xd - t). 

Properties 2. and 3. ensure self similarity and stationarity of the increments of the process. Among 
the kernels satisfying the above properties, [TS] considered: 
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The second kernel defines a class of processes called Hermite processes, which have been and still are 
the subject of considerable interest [51 [251 E 1221 ES]- For d = 1, one recovers fBm (the only SSSI 
Gaussian process). The process obtained with d = 2 is the called the Rosenblatt process. 

When d > 1, the processes in the general class defined in [191 somewhat difficult to analyse, 
because they are no longer Gaussian. In recent years, the family of Hermite processes, and specially the 
Rosenblatt one, has been studied in depth from the point of views of estimation [5], stochastic calculus 
[29j, distributional properties [30j, wavelet-like decomposition for d = 2 [22], and more. Wavelet-like 
decompositions, in particular, are useful for investigating the local properties and providing synthesis 
algorithms. 

General results on sample paths properties of ergodic self similar processes, such as local or 
pointwise regularity or Hausdorff dimensions were obtained in [26]. One such result is that, if the 
process has stationary increments, then for all t, the pointwise Holder exponent is equal to H almost 
surely. In the case of fBm, a lot more is known, both in the one and multidimensional cases [2H I51I27]. 
However, in the non Gaussian case d > 1, there is still room for improvements. For instance, precise 
uniform results are not available. Of course, by stationarity and self-similarity, one has: 



By Kolmogorov's lemma, X has a modification whose sample paths are Holder-continuous of all 
exponents smaller than H (we will always work with such a version in the following), but obtaining 
the exact uniform pointwise Holder exponent is more difficult. 

The aim of this article is to study a class of SSSI processes obtained by considering a particular 
kernel h satisfying conditions 1., 2., 3. above. The interest of this class is that it allows one to obtain a 
wavelet-type expansion for all values of d (and not only for d = 1 and d = 2 as is the case for Hermite 
processes). This expansion permits in turn to deduce sharp local regularity results on the sample 
paths. Our class of processes is defined by the following multiple Wiener-Ito integral representation: 



where t G [0, 1], t* = {t, ...,t) and a = H - 1 + ^. When d = l, this is just fBm. 

We note that, although the class of processes above does not occur as limit processes as is the 
case of Hermite processes, the structure of their kernel is justified by the following formal analogy: 
when passing from one to several dimensions, fBm may be generalized in two ways, leading on the 
one hand to fractional Brownian sheet: 



Hermite processes are the formal analogues of fractional Brownian sheet, while our class would 
correspond to Levy fractional Brownian field. 

The remaining of this work is organized as follows. In the next section, we study basic properties 
of our kernel and define precisely our process. In Section [3l we prove an almost sure wavelet-like 
expansion converging uniformly on every compact. Section HI studies the pointwise and local Holder 
regularity of sample paths and the behaviour at infinity of the process. Finally, we provide in Section 
[5] general results on the Hausdorff dimension of the range and graphs of multidimensional anisotropic 
SSSI processes defined by multiple Wiener integrals. 



Vp > 1, E[\Xt - X,\P] = \t- sr^E[(Xi)P]. 




(1) 




and on the other hand to Levy fractional Brownian field: 
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2 Definitions of the kernel and process 

Let H G (5,1) and a := H + ^ — 1. Consider the Riesz kernel: 



1 jllxll^-'^ 0-^7^0,2,4,6 
^ 1 IMir'^ln^ a-d = 0,2,4,6 



with: 



2"7r2 



a r(f; 



r(^) 



1 



d—a ■•^ 

d-2 d „ , , a — (i , 



a = -2k 



-1)— 7r22"-^(^— )!r(f) a = d + 2k 



where A; G N and T{a) = x" ^dx. This kernel occurs in the definition of the Riesz potential 
which will be used below. Consider the following kernel: 

ht{xi,...,Xd) = ka+i{t* - -x.) - ka+i(yi) 

where t* = {t, ...,t). Let us prove that hf verifies Properties 1., 2., 3. Let t G (0; 1]. In order to show 
that hf is square-integrable, we need to study its behaviour close to and at infinity. Around 0, we 
have, for some > depending only on d: 



|x||2^-'^dx<C^ / r^^-V-^dr < 00 



'B(0,e) ^0 

since 2H — 1 G (0, 1). To analyse the behaviour at infinity, let us first write: 



||t*-x| 



X 



X 



X 2 
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Since 



ii-i^i<"^:,","ii^<^+i 



X 2 



and t* G [0, 1]'^, one has, for ||x|| large enough: 



||t*||2 ||t*-x||2 ||t*||2 



X 2 



X 2 



For d > 2, - f < and 



x,,-ff- 



2 



< 



(1 _ - 1 

X2 



X 2 



||a;||2^+oo 2 ||x||2 



Thus, at infinity: 



Now, 



lit* 



, (^-^)l|t1|2||x||2 



I 



\2H-d-2 



dx. 



PC 



r^H-d-2+d-i^^ < 00 



/]Rd\B(0,e) 

since 2H — 3 G (—2,-1). The kernel hf is clearly symmetric and simple calculations show that 
Properties 2. and 3. are also verified. Thus the process: 

VtG[0,l] X^ = Ia{hf), 
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where denotes the Wiener integral iterated d times, is an SSSI process by [19]. 
Self similarity and increments stationarity entail that 

V(s,t) G [0, E[X^Xf] = ^E[{Xff] + |s|2^ -\t- , 

i.e. has same covariance structure as fBm. In particular, it displays long range dependence. Note 
that, as mentioned above, Kolmogorov's lemma implies that the process has a modification whose 
sample paths are almost surely continuous. We will always work with this modification in the sequel. 
Using a classical patching argument, X°' can then be defined on the whole of M+. 



3 Wavelet decomposition 

Wavelet decompositions of stochastic processes are a powerful tool for the study of sample paths 
properties such as regularity. They have been used in particular in [4J and [2] respectively for the 
cases of fractional Brownian sheet and linear fractional stable sheet. Following the method in these 
articles, we shall obtain an almost sure wavelet-type expansion of {-^f}- 

In that view, we use a specific wavelet basis of Li^iW^) (the reader is referred to \VJ\ for an 
introduction to wavelet theory). For d = 1, we consider the Lemarie-Meyer wavelet basis, which 
verifies the following property: 

• the scaling function (p and the mother wavelet -0 are in the Schwartz space, ^(M). 

• supp{ii)) C G M : 1^1 < ^} and supp{ij) C G K : ^ < |^| < 

• {22^{2^x - k) : {j,k) G Z x Z} is an orthonormal basis of L (M.). 

• V/3 G N Jj^x^'ip{x)dx = 0. 

The last property implies that the mother wavelet is in the Lizorkin Space, S'o(M). We use the 
canonical construction of a multidimensional multiresolution analysis from a monodimensional one 
using tensor products. Let E denote the set {0, l}"^ \ (0, 0). For all e G E, we note ip^^^^ the function 
defined by: 

where = 4> and ^^^^ = V'- Then {2V^/;W(2Jx - k) : (j,k,e) e Z x Z'^ x E} is an orthonormal 
wavelet basis of L^(M'^). By construction, for all e £ E, ijj^^^ is in the following Lizorkin type space: 

5o(M'^) = {ij£ S(R'^) : V/3 G N'^ / x^ij{x)dx = 0}. 

This space is invariant by the Riesz potential which is defined by |23) : 

rV'(x)= / A:„(x-y)V^(y)dy. 

For ip in S'o(M'^), the following inversion formula holds: 

We are now ready to give our wavelet- type expansion for the process X^. 
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Proposition 1 Almost surely, 

yt e [0,1] = E E [^"+'(^^^^)(2^'t* - fc) /d(V'g). 

iez fcezd eeE 

In order to show this proposition, we need the following lemma: 

Lemma 2 Let {fn} be an orthonormal basis of iF'iW^). There exists an event Vl* of probability 1 and 
a strictly positive random variable, Cd{u}), which has all its moments finite such that: 

Vu; G 17* Vn G N \Wn)\ < (7rf(a;)(log(e + n))i 

Proof: The proof follows the lines of the one of Lemma 1 in [3]. Let a„ = (log(e + n))2 . Theorem 
6.7 of [13] entails that: 



3cd > : V6 > 2, Vn € N, F (\Idifn)\ > banVdJ.Wfnh) < e-''^^'"'-^' 



(e + n)'=d&3 



For Cfibd. > 1, ^ — ^ < oo, and the Borel-Cantelli lemma implies that: 

3n*, F{n*) = 1 : Vw G f)*, G N : Vn > n{uj), 

|/d(/n)M| < 6(log(e + n))iVd!||A||2 < 6(log(e + n))i^/d!. 



Set: 

[min{n G N : /c > n |/d(/fc)| < 6(log(e + A:))i\/d!} w G 17* 

1 +00 u ^ 0*, 



r(u;) 
and define: 



sup i\Wk)\) oj^n* 

C(a;) = { o<k<T{uj) 

+00 oj ^Q*. 

We compute the moments of any order p > of this random variable. Following the proof of Lemma 
1 of 0, one gets: 



Ei\cn = Y.^i sup mfk)\nx{rH=n}) 

+00 n 
+00 n 

< E{\Wo)n + Y,J2^mfk)\'n'mWn-i)\ > KMe + n - l))f 



'{|Jd(/n-i)|>6(log(e+n-l))^v^}^ 
n=l k=0 



n=l fc=0 

Using Theorems 5.10 and 6.7 of [1^, one deduces 



E{\c\n < maifoW) + cpE — 

(e + n-l) 



n=l 



2 ' 



i.e., for > 2, 

Vp>0 E(|C|P)<+oo. 
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From this, one deduces that there exists a strictly positive random variable Cdiui), whose every 
moments are finite and such that: 

a.s. VnEN < QH(log(e + n))f □ 

Remark: The main argument used here is of hypercontractivity type. Indeed, the hypercontrac- 
tivity of the Ornstein-Uhlenbeck semigroup allows to show that all the norms are equivalent on any 
Wiener chaos spaces. Since, for all n, Id{fn) is in the Wiener chaos of order d, we can apply that re- 
sult to this quantity. For further details, we refer the reader to the chapter 5 of [13] or chapter 1 of |20j . 

An indexing argument shows that the following inequality holds: 

a.s. V(j,k,e) G ZxZ^xij; \Id{i^^l)\ ^ CiA^){\og{e+\j\ + \Mi))-2 < (5i,d(c^)(log(e+|j|))i(log(e+||k||i))i, 
where we have used the two trivial facts: 

• the function x — )• X2 is increasing on 

• V(x, y) £ M+ log(e + x + y) < log(e + x) log(e + y) 

Let us now move to the proof of Proposition [TJ It is similar to that of the wavelet-type expansion 
of fBm in [IS]. 

Proof: Let t be in [0,1]. Since hf is in L^(M'^), we can expand it into the wavelet basis {V'j'^k}' 
Using Ito isometry we have, in (L^): 

Let us show that the right hand side converges uniformly in t on W , i.e. that: 

y^oen* E E Ell < ^^^S > WoolWfl)] < +00 

For this purpose, write: 

< hf,ipfl >= ^ / [A;,+i(t* - x) - A:,+i(-x)]v]^i(x)dx. 
The change of variables — k = u yields: 

< hf,M >= 2-^'^ . \ ^. [ [K+i{2h* - k - u) - K+ii-k - u)]V'(^)(u)du. 
Split the integral into two terms to get: 

j"+i(^(e))(2Jt* - k) - I''+\4'^){-k) . 
We need to distinguish between j < —1 and j > 0. Since 

g SoiM.'^), the function / defined 

on M by: 

/(n)=r+i(V^(^))(u*-k) 
is C°°(M). Thus, by the mean value theorem: 

3v G [0,2-''t] : /"+i('0('))(2Jt* - k) - /°+i(V'('))(-k) = 2H^{v). 



<h?,i;fl>=2-^ 
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By the chain rule, 



i=l 



Using the rapid decrease of one gets: 



Vp>0 \^{v)\< ^"''^ 



'du^ " ~ (1 + ||v* - k||2)P 
For j < -1, u G [^,2H] C [0, \]. Thus, for ||k||2 large enough, the following estimate holds: 

By Lemma [21 

V^^ G f]* V(i,k,e) € Z X Z'^ X \h{^%){^)\ < (5(a;)(log(e + U|))i(log(e + ||k||i))i. 

Thus, 

I < > ||/,(^ t)H| < (log(e + |j|))Mlog(e + ||k||i))2. 

(I + II2 -k||2F 

Since H < 1, almost sure uniform convergence for the "low frequency" part j < for p sufficiently 
large does hold. The high frequency part j > requires a bit more work. The inversion formula for 
the Riesz potential yields: 

j"+i(^W)(2Jt*-k) = ^ [ e*<2^**-'^'«>^%de 

Making p integrations by parts in the direction i where \\2^t* — k||2 < Cd\2^ — ki\, one gets: 

c 

' ^' ~ (1 + 2n* -k 2>' 



1 

+ 



Thus, for j > 0: 

<n,,^.^j,> I SS.rf^ [(l + ||2it*-k||2)p ' (l + ||k||2)p 

In order to conclude, we need to estimate the following quantity: 

^ (log(e + ||k||i))i 
(l + ||2Jt* -k||2)P' 

Denote [2-'tJ* = {[2H\, [2H\) G Z'^. A change of variable in the previous sum yields: 

^ (log(e + ||k||i))i ^ (log(e + ||k+ L2^tJ'-||i))i 
^ (1 + ||2it* - k||2)P ^ (1 + ||2it* - |2n|* - k||2)P' 

By the triangular inequality, 

log(e + ||k+ L2^'tJ*||i) < log(e + ||k||i + ||L2^tJ*||i) < C'log(e + ||k||i)log(e + ||L2^tJ*||0 
Consequently, one obtains: 

y < c'(iog(e+iiL2^-tj*iio)i sup ( y m^+\m^))^ 

^^ji + \\2n*-kUp- J iii^\^j^P,(^Z^^ (i + l|,_k||,). 

Thus, for j > 0: 

E I < ^"'^S > < C;;,(a;)2-^-^(log(e + d2^-t))i(log(e + 
which concludes the proof .□ 
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4 Uniform and Local Regularity of the sample paths 

From Proposition [U one sees that the process may be written: 

where X^'" is the low frequency part and X^''^ the high frequency part, namely: 

In order to study the regularity of the stochastic process, we will show that the low frequency part is 
C°°(R) and then focus on the regularity of the high frequency part. 

Proposition 3 Almost surely, the process {x''"} is C°°(M). 

Proof: Let gj]^ be defined by: 

VtG[0,l] g,^^,{t) = r+\^l,^^){2H* -\,). 
We need to show that the following random series converges uniformly on [0, 1] almost surely: 

^m„ 

E EE2-"^(^)Wti), rn^W. 
The chain rule formula entails that: 

Since G 5*0 (K'^), one obtains the following estimate: 

I " yi.k/,\| ^ r,jm '^p,d,m 



VP>0 \-j^{t)\<2^' 



dt"^ ^ ~ (1 + \ \2n* - k\\2)P' 
Moreover, since t G [0, 1] and j < -1, 2h* e [0, ^f. Thus, for ||k||2 sufficiently large, one has: 



dt^ - (i + ||i*_k||2)p' 

Finally, since H < 1 < m and using Lemma [21 one obtains the almost sure uniform convergence of 
the random series. □ 

The next two results are analogues of Theorems 1 and 2 of [5]. The first one gives a modulus of 
continuity for the process X" while the second one deals with the asymptotic behaviour of the pro- 
cess as t — )• +00. Their proofs are also similar. 

Proposition 4 There exists a strictly positive random variable Ad of finite moments of any order 
and a constant, bd > 1, such that: 

v.en- sup \IlM^m:^ < A,i.). 

(s,t)e[o,i] \t - s\^(\og{bd + \t- 
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This result should be compared with Remark 2.3 p. 166 and the sentence that follows its proof on 
p. 167 of [26], which provide a related (but weaker) result, namely that the local Holder exponent of 
X" (as well as the ones of all Hermite processes) is equal to H: the above proposition entails that 
the local Holder exponent of X" is not smaller than and Theorem [10] below will ensure that it is 
not greater than H. See also Lemma 5.9 of |25] . which provides a similar (though slightly weaker) 
result for Hermite processes. 
Proof: 

Since t — t- x''" is in C°°(IR) almost surely, we need only to deal with the high frequency compo- 
nent X^'°^ . Let a; G ri* and (s,t) be in [0, 1]. There exists a positive integer jo such that: 

2-^°-^ <\t-s\< 2-^0. 

Write: 

= (/) + (//), 

where 

30 

(/) = E E [^"+'(V'^^^)(2^'t*-k)-r+i(v(^))(2^s*-k)] w];i) 

and 

(^^)= E EE2"'''[^""''(V'^'^)(2''t*-k)-/"+'(V'(^))(2^s*-k)]/,(^jJ). 
i>io+i keZ'* £££ 

Let us consider the first sum. Using Lemma [2] and following the lines of the proof of Theorem 1 of 
[1], there exists a deterministic constant hi^d > 1 such that: 

io 

1(1)1 < A,^d{u:)\t - ,| J]2(i-^)^-(&i,, + j)f < A^A^)\t - s|2(^«+i)(i-^)(6i,, + jo)i 
i=o 

Similarly, for the second sum, there exists a deterministic constant 62, d > 1 such that: 

00 

|(I/)| < A2,d(t^) 2-^'^(62,d+j)^ <^'2,,(^)2-^«^(62,d + io)^. 

i=io+i 

Putting the results altogether and using the definition of jo yields the desired conclusion. □ 



Proposition 5 There exists a strictly positive random variable Bd of finite moments of any order 
and a constant > 3, such that: 

Vw e n* sup l^lMl _ < Bdico). 

tm+{i + \t\)H{loglogicd + \t\))-2 

Proof: The proof is quite similar to the previous one except that the behaviour at infinity of the 
process will be governed by the low frequency part. Let t be a real strictly greater than 1 and let ji 
denote the integer such that: 

2-'i < t < 2-'i+^ 

Let uj G Q*. We consider the following quantity that we split into three parts: 
E E E^"'"" [l''^\r)i2H* - k) - /-+HV'^)(-k)] (log(e + \j\ + ||k||))f = (/) + (//) + (///) 
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where, 

-(ii+i) 

(^)= E E ^2-^^[l'^+\i,^)i2H*-k)-I-+\r)i-^)] (log(e + |i| + ||k||))i, 

(//) =f^Y.H 2-^-^ - k) - r+^(V^)(-k)] (log(e + \j\ + ||k||))i, 

j=-Ji kezrf eeE 

oo 

(m) = ^ E E 2"''' [I^^'mi^'t* - k) - J"+\V^)(-k)] (log(e + lil + ||k||))i. 
Let us deal first with the high frequency part. We have: 

|(m)| </i3(t) + /i3(o) 

where is the function defined by: 

oo 

^3(*) = E E E2"'''l^"+'(V'^)(2't* - k)|(log(e + IjI + ||k||))i. 
j=o kezd eeE 

Reasoning as above, one gets: 

+00 ^ 
hsit) < C3^2-^-^ (log(e + J + ||2^t*||))^ . 

j=0 

Using that: 

Va, 6 > log(e + a + b) < log(e + a) + log(e + b), 
and, for d > 2, by a convexity argument: 

Va, 6>0 {a + b)^ < Cd{a^ +b^), 

one gets: 

+00 +00 ^ 

hsit) < Csd Yl 2"''' (log(e + j))-^ + Cm ^ 2-^" (log(e + 2^||t*||)) ^ . 

Finahy, since log(e + 2^\\t*\\) < log(e + 2-') + log(e + ||t*||): 

hs{t)<C',a(.log{e+\\t*\\))t 
Consider now the second term. One has: 

\{n)\<h2{t) + h2io), 

where /i2 is the function defined by: 

-1 

h2{t)= E ^^2-^-^|r+HV'(^^)(2^t*-k)|(log(e + |j| + ||k||))i. 

j=-ji keZd eeE 



It is clear that: 



h2 it) < C2 ^ 2^-^ (log(e + J + 2-^ 1 1 1* 1 1 )) ^ 
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Using the definition of ji and elementary computations, one obtains: 

/l2(t) < C2,d{l + \t\)^ (loglog(crf + \t\))^ . 

Similarly for /i2(0), one has: 

h2{0) < + \t\f (loglog(c, + \t\))-2 

For the last term, we apply the mean value theorem to the function / defined above to get: 

j=-oo keZdeeE ^ " "'^ 

and the result follows. □ 



Remark 6 In order to study the local regularity of the process, we will only need the behaviour at 
infinity of the high frequency part. 

We now compute the uniform almost sure pointwise Holder exponent of X°^. Since the paths of 
the process are nowhere differentiable, it is by definition equal to: 

Xa va I 

I t-\-o — t I 

7X"(i) = sup{7 > : limsup-^ p-j < +00}. 

Proposition m entails that, almost surely: 

VtG(0,l) ^x4t)>H. 

We will show the converse inequality. In that view, we introduce a stochastic field X" that is closely 
related to Xf and defined as follows: 

v.. G n*, vt G R'^ x^uj) = E E E^"''' [^"+'(V'^^^)(2^t - k) - r+nv'^^^)(-k)] hi^jM). 

It is easily checked that X" is well defined and that its sample paths are continuous. Moreover, 

a.s., Vt G |X^| < (7d(a;)(log(e+ ||t||))i . 
Define the following operator on 5*0 (M'^): 

VxGM'^ D-+i(^)(x) = / e^<"^€>||eir+^V'm. 
The Fourier transform of the dilated and translated version of this function is equal to: 

where J-(Z)"+i('i/'))(6 = IICIT+^V'IC)- 
Lemma 7 Almost surely, 

V(j, k, e)N xZ'^xE Idii'fl) = 2^'^^+'^^ [ X'^D"+\^P'){2H- k)dt. 
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Proof: For w G 0* and (j, k, e) G N x Z'^ x we compute the following quantity: 
A:= [ Xf^ ^^D''+\^l;^'^){2Hl-h,...,2Hd-kd)dh...dtd. 

jRd. 

Since the function D'^^^{ilj^''^) is in 5*0 (M*^), the asymptotic growth properties of the stochastic field 
{Xf } allows one to invert the sum and the integral, so that: 

A = V 2-^"^ /" /"+i(V'(^'))(2^"t-k')I)"+^('(/'^'))(2^t-k)dt /d(^^fb 
=: J2 ^i'>k',.' Wfl')- 



Let us compute the inner products A^, ]^^^,. By Fourier isometry, 



Using the inversion formula and the definition of the operator Z)"^^, one gets: 

r)-j{H+d) /-,,./ . . 



Since: 
one finally obtains, 

But {V'j'^k} orthonormal basis of L^(M^), so that: 

f x^L»"+^(V'')(2^'t - k)dt = 2-^'^^+'^hd{%bfl). □ 

We will show that the pointwise Holder exponent of the stochastic field X" is everywhere equal to H 
almost surely. As a first step, we will exhibit an event of probability one such that, for any to G [0, 1]*^, 
there exists a random variable Id{i^^j\^ which is both "large" enough and whose index is "close" to 
to- In that view, we will use the following result from [6]: 

Lemma 8 |3, Theorem 8] There exists an absolute constant c > such that for all polynomial 
Q : M" — )• M o/ degree at most d, all 1 < q < oo, all log-concave probability measures fx on M" and all 
a > 0, 



\Q{x)\d iji{dx) J /x({a;GM" \Q{x)\ < a}) < cqa^ . 

Lemma 9 Let 1 < /3 < 2. There is an event il** of probability one such that for all oj £ 17**, for all 
e in E, for all to G [0, 1]°', and for j large enough, there is a kj G Z"' such that: 

\to-2-^kj\ < bj^2-^ 

where b > is a constant independent of co and 

\Wj,k,)\ > V 

where rj > is a constant independent ofuj. 
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Proof: Let j be a positive integer. For simplicity, we assume e = (1, 1). The general case follows 
in a similar way. We consider the dyadic points of order j of (0, 1)°' i.e. {^.2^^ : k G {1, 2^ — 1}'^}. 
Set Aj = {1, 2^ - ly. Let us build a covering of A^. In that view, define the following sets: 

VqGQ,, Dj = {r/lq + r:rG{0,...,r/l -in. 

Clearly: 

A, C U Df. 

Consider the following event: 

Where diag{Dj) denotes the diagonal of Dj. We wish to estimate its probability, and in particular its 
dependence on j. For this purpose, we first perform a study of the random variables -^d(V'jk)- Since 
e = (1, 1), one has: 

and 

Let us reorder the previous symmetric tensor product in the following way: 

where I G 1, d, 71 + ... +7; = d and, for all p, m in {1, I}, kp ^ km. Using the formula expressing 
Wiener-Ito integrals in terms of Hermite polynomials, one obtains: 




The standard normal random variables ji'^ipj^kp{u)dBu are independent. As a consequence, each 
Idii^j]^ is a polynomial of degree at most d evaluated at independent standard normal random 
variables. Moreover, for each k,k' G diag{Dj), Idii^jk) and /d(V'jk') mutually independent for 
k ^ k' . Therefore: 

Lemma [8] with q = 2d yields 
Using the trivial fact: 

^^Il^j,kll2>l, 

one obtains the estimate 

n{\wik)\ < < 

which entails that 

2^ , 1 fil^l 

P(A,)<(L-^J+2)'^(-)V. 
13 



Since /? > 1, Yli=i ^i^j) < °o- By the Borel-Cantelli lemma, the fohowing event has probabihty one: 



oo oo 



j=ij=j 

where denotes the complementary of the event Aj. Let lo £ Q** and to G [0, l]'^. There exists 
kj G Aj such that |to — kj2~-' | < Moreover, since UgeQ ^'j ^ covering of Aj there exists qj 

such that kj G D^^. But for j large enough, there also exists kj G such that \Id{i^ji.)\ ^ 
Write: 

|to - 2"%j| < |to - kj2"^'| + 2"-''|kj - icjl < 6/2"^ 
which yields the desired conclusion. □ 

We are now ready to prove that the pointwise Holder exponent of the stohastic field X" is equal 
to H almost surely: 

Theorem 10 Almost surely, 

ViG(0,l)'=' ^x-{t) = H. 
Proof: Assume that there exists ojq G il**, to G (0, l)'' and eo > such that: 

7x«(to)(^o) >H + eo. 

By definition of the pointwise Holder exponent, for any t such that ||t — to||2 is sufficiently small, the 
following inequality holds: 

|X?(a;o)-X-(..o)| <CMI|t-to||f+^«, (2) 

where C{ujq) > is a finite constant. By continuity and growth at infinity of X°, this inequality 
holds for any t G W^. For any (j, k, e) G N x Z'' x 

|I,(V]k)(a;o)| = 2^^(^+'^)| / X^Z)"+i(^^)(2^t - k)dt|. 

Using the fact that G Sq{W^) and inequality ([2]), one gets: 

< 2^-(^+'^) / \X^-X^J\D-+\^P^)i2H-k)\dt 

< C7(a;o)2^'(-^+'^) / ||t - to||f +'''|L'"+^(V'')(2^'t - k)|dt. 
The rapid decrease of D'^^^{Tp'^) and the the change of variable x = 2-'t — k yield 

|/d(^-,k)K)| < C'(a;o)2-'-^ / ||2--^x + k2-^--to||f+^° ^ ^,1 ,, . dx. 

Jk^ (i + ||x||2r 

Since (|t| + |s|)^+^o < C"{\t\^+'° + |s|^+^o), one obtains 

k)K)| < C'"(a;o)2-^'^«(l + ||k - 2^to||f +^»). 
This inequality is true for any j and k. Applying it with the kj of Lemma O yields 

\Wl^){u:o)\ < C""K)2-^'^°(1 + (6/)^+^"). 

This right hand side in the above inequality tends to as j tends to +oo. But, by Lemma O 
{Iditpj ]^.){ijJo)\ > V- Then, ad absurdum, 

a.s, Vt G (0,1)'^, 7x-(t) = H. □ 
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Corollary 11 Almost surely, 

Vt G (0,1), jxAt) = H. 

This result should be compared with Proposition 2.2 p. 163 and Remark 2.1 p. 164 of which entail 
that, for all t, almost surely, the pointwise Holder exponent of X" (as well as the ones of all Hermite 
processes) are equal to H. 

5 HausdorfF Dimensions of SSSI processes represented by multiple 
Wiener-Ito integrals 

In this section, we provide estimates for the Hausdorff dimension of the range and the graph of 
multidimensional processes represented by multiple Wiener-Ito integrals of order d whose kernel verify 
conditions (1) — (3) with self-similar exponent H. More precisely, let d > 1 and A'^ > 1 be two integers. 
Let {Y/^} be the process defined by: 

Yf = -f{H,d)I,{h^) 

where 7(-ff, d) is a normalizing positive constant such that EdYj^p] = 1. 

Let ^ < Hi < ... < Hn < 1- Let {Y^} be the multidimensional process defined by: 

{Y^} = {{Y,"\...,Y/'-):t€R+} 

where the coordinates are independent copies of the process Y/^ . 

Let E C M+. It turns out that the methods developed in [31J apply with few modifications 
for estimating the Hausdorff dimension of the range REiY^) := {Y^{t), t e E} and the graph 
GteO^^) := {{t,Y^ {t)), t ^ E} oi Y^ over E. We shall make use of the following lemma: 

Lemma 12 There exists a stricly positive constant, C , depending only on d, such that: 

Vj e {l,...,iV} ¥{\y"'\<x) <Cx^. 

Proof: We shall use Lemma[8]in order to prove this inequality. For all j G {1, A^}, ^{Hj, d)h-^ ^ G 
L^(M'^). If {fn -nG N} is an orthonormal basis of L^(M), then {fi^ ... fi^ : (ii, id) S N'^} is an 
orthonormal basis of L^(M'^). Thus, one can approximate 'y{Hj,d)hf^ by: 

n 

^f""= Yl <l{Hj,d)hf';fi,(S)...(g)fi^>fi,(g)...(S)fi^. 

ii,...,id=l 

By linearity of the multiple Wiener-Ito integral, 

n 

/d(/if"") = <l{Hj,d)h^'-fi^®...®fi^>Id{fn®...®fi^). 

il,...,id=l 

The Wiener-Ito theorem entails that /^(/n ••• fid) is a polynomial of degree at most d evaluated 
at independent standard normal random variables. Therefore, 

/d(/if^'") = P„,d( / fi{u)dB,{u),..., [ U{u)dBj{u)) 

JR JR 

where P^.d is a polynomial of degree at most d. Applying Lemma [8] to this quantity with q = 2d and 
H the Gaussian probability measure on M" yields that there exists a constant C{d) > such that: 

VneN Vx>0 P(|/d(/if"")| < x) < C(d)x5||/if^'"||2^. 
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By the continuity of the multiple Wiener-Ito integral, Idih^ ^' ) — )• Id{^{Hj, d)li^^) in L'^(i7). Consider 
a subsequence such that Id{hi ^' ) converges almost surely to Id{j{Hj, d)h^ ^). Using Fatou's lemma 
along this subsequence, one obtains: 



P(Id(/ifO <x) < liminf P(|Id(/if^'")| < x) 
< C{d)x^ liminf H/if^'^IU 

s>+oo 



H 



< C{d)xd \\-i{Hj,d)h{ 

< C{d)x^{^)-^ 

Wdv 



J I I d 
2 



< K{d)xd 

since, we have ^[\y"'\'^] = -i{Hj,dfd\\\h^'\\l = 1. □ 



We are now ready to prove the following theorem which is similar to Theorem 3.3 of [3T 
Theorem 13 Almost surely, 



r- p I Y!j = liHk-Hj) 

dimuREi^") > min ( N; ^ , k = 1, ...,N \ , 

Hi. 



dimnGrE{Y^) > mm \ ^ ,k = 1, ...,N,dimnE + '^^ ~ 



Hh. 



1=1 



d 



Proof: Assume that dimy^E > (otherwise there is nothing to prove) and let k be such that: 



V^fc— 1 TT TT 



d 



d 



Let us compute the quantity: 



dimruE + ^^=1*-^^' 
A = min I N; ^ — ^ ,j = l,.:,N 



Let j < k. Let: 



A := 



dim-uE + ■^') dimnE + ^0 



Then, 



A 



1 



HjHk 
1 

HjHk 



dim-}iE{Hj — H^) + 



Yli=i(.Hk - Hi)Hj Y^{=i{Hj - Hi)Hk 



d d 



d 



d 



Since: 



k k-l k-l 

(Hk — Hi)Hj = ^ (i/fc — Hi)Hj < ^ (fffe - Hj)Hi, 

i=j+l i=j+l *=i+l 
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one gets: 



Hk 



HjHk 



d 



which is strictly less than by definition of k. One proves in a similar way that this inequality also 
holds for j > k. Thus, 



A 



From this one deduces that < A < ^. 



Hk 



We shall use a classical potential theoretic method to obtain a lower bound on the Hausdorff 
dimension. More precisely, in the case of the range of {Yj^}, we study the convergence of the following 
quantity: 

dfiYH{x)dnYH{y) 



where 7 is a strictly positive constant and is the image by of a positive measure, on E, 
both of them suitably chosen. By a change of variable, the previous quantity is equal to: 

dii{t)d^ji{s) 



By Fubini's theorem, this reduces to the convergence of the following integral: 

/ E[|Yf -Yf|-^](i/i(t)dMs). 

JExE 

Fix 7 e A) and set 



TV 



F{x) = JYA\Yi\ < 



One computes: 



E[|Yf-Yf|-T] = [ [{xi\t-s\^'f + ... + {xN\t-s\^^f]-^dF{x) 

= \t-s\-^'^ [ dF2{x2)...dFN{xN) [ [xl + ... + {xN\t-s\"''-"'f]-^dFi{xi). 

For all a > and all /3 > 0: 

J := [ [xl + a^j-^dFiixi) = /3 [ xiFi{xi)[xl + a'^]-2-'^dxi. 



When f3 > ^, Lemma [T2] yields: 



+00 



J < f3C{d) I [xl + a^j-^^^xl^'^dxi 



< C{d) 



JO J a 



< C'{(3,d)a-^+^. 
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When ^ < i, 

"1 01,1 /•+00 



and thus 



Jr+ Jq Jl 

[ [x2 + a2]-ldFi(xi)<C7"(/3,d). 



Since7 G (^^;-^) we see that, upon integrating iteratively with respect to ciFi(xi), dF2(x2), (ii^A;(xA;), 
the integral 

[ [{xi\t - s\"'f + ... + {xN\t - s\^^f]-2dF{x) 
is not larger than a constant times: 

\t - s|-^^i7-(^^2-Hi)(7-^)-...-(Hfe-H,_i)(7-^) ^ 1^ _ 5|-^^fe7+E-=i 

By Frostman's Lemma, the condition -f^fc7 — X]j=i — ~ ^ dimy^E entails that there exists a positive 
measure fi on E such that 

^(ds)/u((it) 



< +00. 



Therefore, 



E 



fi{ds) fi{dt) 



< +00 



I.e. 



sxi? |Yf -Y|^|7 
dimnREij") > 7 a.s. 



for aU 7 G (^,A). 



The proof of the lower bound for the Hausdorff dimension of the graph of over E follows from 
the one of Theorem 3.3 in [31] with the same modifications as above. We omit the details. □ 
Remarks: 

• The proof relies essentially on Lemma [T2l The lower bound holds for multidimensional Hermite 
processes of any order d as well as for multidimensional versions of X". 

• By studying further the law of , one would obtain a better lower estimate. In particular, it 
is well-known that the density of the Rosenblatt distribution is continuous and bounded ([1], 
[7]). As a consequence, P(|y^ ^ | < a;) < K{d)x, which yields the following lower bound: 

/ dimnE + y^^-AHk- HA \ 
dimnREi^ ) > min ( N; -,k = 1, N j a.s. 

( dim-uE + y^^-AHk- HM \ 
dimnGrE{Y") > min ( k = 1, ...,N, dimnE + ^(1 - Hi) j a.s. 

One thus retrieves here the result of Theorem 2.3 of |31j . 

Upper bounds for the Hausdorff dimensions of the range and the graph are obtained by a direct 
application of Theorem 3.1 in [31], since we deal with elements of Wiener's chaoses, which thus have 
finite moments of all order: 
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Theorem 14 !131\ Thoerem 3.1] Almost surely, 

(iim^/?ij(Y^) < mill I iV; ^"^^ ^,fc=l,..,iVj , 

dimnGrEiY^) < min ( k = 1, N; dimnE + ^(1 - Hi) 

Remarks: 

• For the A^-dimensional anisotropic Rosenblatt process, we retrieve the known resuhs: 

R fwH, ■ (m dimnE + ZU(Hk-H,) \ 

dimnRE\y ) = mm A*; ,k = 1, ...,N a.s., 

\ Hk J 

( dim-HE + y^''_AHk- HA ' 
dimnGrE{Y^) = min ( k = 1, N; dimnE + ^(1 - H,) 



a.s. 



• In the isotropic case, the first equahty was proven in |25j . 
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